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Abstract 

Recent Couette-cell shear experiments of carbopol gels have revealed the formation of a tran¬ 
sient shear band before reaching the steady state, which is characterized by homogeneous flow. 
This shear band is observed in the small-gap limit where the shear stress is spatially uniform. An 
effective-temperature model of the transient shear banding and solid-fluid transition is developed 
for the small-gap limit. The small-gap model demonstrates the ability of a continuum-constitutive 
law that is based solely on microstructural rearrangements of the gel to account for this tran¬ 
sient behavior, and identifies that it proceeds via two distinct processes. A shear band nucleates 
and gradually broadens via disordering at the interface of the band. Simultaneously, spatially 
homogeneous fluidization is induced outside of the shear band where the disorder of the gel grows 
uniformly. Experimental data are used to determine the physical parameters of the theory, and 
direct, quantitative comparison is made to measurements of the structural evolution of the gel, its 
fluidization time, and its mechanical response under plastic flow. 
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I. INTRODUCTION 


Yield stress fluids (YSFs) are ubiquitous in everyday life, and their special properties have 
merited intense research |lR9]. Examples include gels, clay suspensions, foams, concentrated 
emulsions, and colloids. These seemingly distinct substances exhibit a similar mechanical 
response when subjected to shear deformation: Below a critical (yield) stress these materials 
remain elastic, behaving as solids, but above this critical stress they are able to deform and 
flow as viscous liquids. This characteristic ability makes them extremely sought-after for 
many applications 10]. In the presence of an applied shear protocol a YSF can exhibit a 


distinctive mode of deformation known as shear banding , where a region of highly localized 
strain that is far greater than that of the surrounding material appears. The region, or band, 
is often long or continuous in the direction of shear and of some finite width perpendicular 
to the shear. This form of shear localization has been shown to be a common occurrence in 
the rheology of complex fluids as well as glassy materials in general fll |. 

YSFs are often categorized by the degree of thixotropy they present. Thixotropic YSFs 
have memory and aging effects that lead to a history dependence, marked e.g. by flow curves 
that are significantly different when the applied shear rate is ramped up compared to the 
ramp down. In contrast, nonthixotropic or “simple” YSFs have no apparent history depen¬ 
dence over the timescales of observation. They also tend to have predominantly repulsive 


interactions between their constituent mesoscopic substructures 


12]. More important is the 


fact that thixotropic YSFs, regardless of the particular material, easily form shear bands. 
Some known as viscosity-bifurcating YSFs, display shear banding in the steady state, a p he- 


13 


and wormlike micelles 


161 


nomenon which as been observed in a number of complex fluids including polymers 

17 -J 20 ]. Experiments of simple YSFs on the contrary, report conflict 


ing accounts of the formation of shear bands: I 11 some instances shear bands are observed 
and yet in others the system (sometimes the same system) remains completely homogeneous 


under shear 


29] . These contradictory findings have prompted investigations of simple 


YSFs with greater attention to the specific material and composition, the shear procedure, 
and the influence of the geometry of the apparatuses. 

Carbopol gel is widely regarded as a quintessential simple YSF and as such has generally 


been be 
system 


ieved to transition from the solid to the liquid state uniformly as a homogeneous 
, [3]. Recent Couette-cell experiments of shear in a carbopol gel have instead 
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revealed the conspicuous formation of transient shear bands before the gel reaches a uniform 


steady state characterized by a linear velocity profile and Herschel-Bulkley rheology 


26]. 


More importantly because of the small width of the particular Couette cell, the shear stress 
measured in these experiments was observed to be nearly uniform across the gap, and the 
stress gradient was reported to not play a significant role in the formation of the transient 
shear band. 

The onset of plastic flow in YSFs and other similar materials has been investigated 
by other theoretical approaches, but direct, quantitative comparison has been lacking with 
regard to shear banding and the pro cess of fluidization. Among these theoretical frameworks 
are soft glassy rheology (SGR) |30|. The SGR model in particular has been successful in 
qualitatively describing the general steady-state power-law behaviors seen in experiments, 
and criteria have been proposed to characterize the onset of shear banding 3l|. It has also 
been able to capture aging. One drawback of SGR is that it assumes the existence of a 
“noise temperature” that controls activation rates of plastic processes. The physical basis 
for this noise temperature is however not clear, and its dynamics have not been derived from 
fundamental or mesoscale principles. 


A minima 


leoretical model of a Newtonian fluid using the Krieger-Dougherty constitu¬ 


tive relation 32] has been used to qualitatively model transient shear banding in a Couette- 


cell geometry. The shortcomings of this approach are that it does not describe a yield stress 
or solid-state elastic response, nor does it quantitatively capture key experimental observa¬ 
tions of shear banding in YSFs, i.e. the correct spatial and temporal evolution of the shear 
band and fluidization times of the system. More importantly, this model attributes transient 
shear banding to the nature of the Couette-cell apparatus itself. Specifically it suggests that 
the gradient in the stress held, which exists simply due to a cylindrical geometry, leads to a 
higher shear rate and unjamming of material near the inner cylinder of the Couette cell. At 
the present time however, it is unclear to us whether the primary mechanism by which such 
transient shear banding occurs is the stress gradient caused by the cylindrical geometry, or 
if the mechanism is an instability in the microstructure of the gel itself. In the case of hard- 
sphere colloidal glasses similar questions have been raised and subsequent analyses have been 
preformed to address the possibility that an intrinsic instability via shear-concentration cou¬ 
pling can lead to shear banding in the absence of gradient-driving nonuniform stress fields 


33 ]. 
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Motivated by the possibility that heterogeneities in the microstructure of the gel may 
indeed be the dominant driver of the observed shear band, we propose an alternative de¬ 
scription of the phenomenon of transient shear banding during plastic flow based on fluc¬ 
tuations in the structure of the gel. We adapt the effective-temperature hypothesis of the 
shear transformation zone (STZ) theory as a rheological model for shear localization in a 
simple YSF. In this description, structural changes described by an effective temperature 


34] account for the observed plastic flow, and the formation of a transient shear band is the 


primary mode of deformation of the carbopol gel. This occurs simultaneously alongside a 
distinct, uniform fluidization of the gel outside the shear band. We then make quantitative 
comparisons of stress-time behavior and velocity profiles with recent small-gap Couette-cell 
experiments of a carbopol gel. 

The STZ theory is a general framework for characterizing plasticity in amorphous materi¬ 
als, and provides a continuum-level mean-held approximation for how based on an assump¬ 
tion of local rearrangements of a material’s structure via the activation of STZs 34h46]. 
STZs are orientational point defects that mediate plastic how by accommodating rearrange¬ 


ment. In this paper we use data from the aforementioned experiments in 


26] to determine 


the physical parameters of the theory and make precise experimental connections. A dis¬ 
tinguishing feature of the STZ theory is that it is based on a specihc model of molecular 
rearrangements, which have been observed directly in numerical simulations and analog ex- 


, |43, |48|. The presen 


periments (43 


the original work of Argon 


49 


formulation of the STZ theory has been extended from 
50| and his proposal of “shear transformations” to explain 


plastic deformation in metallic glasses, as well as from the free-volume and how-defect the 


ories of Turnbull, Cohen, Spaepen, and others (51 


52]]. The STZs themselves have been 


postulated to have internal degrees-of-freedom |43|]. The STZs not only transform from one 
orientation to another; they are created and annihilated during configurational fluctuations 
at a rate that is proportional to the rate of energy-dissipation per STZ [36]. A coarse-grained 
continuum description of the effect of the changing population density and spatial distri¬ 
bution of the STZs is dehned by an effective-temperature held. In this paper we present 
equations of motion for the effective temperature and plastic-strain where the steady state is 
set by a Herschel-Bulkley rheology, and no explicit dependence of the STZ internal variables, 
local strain rate, or how history appears. 

A how rule for the plastic component of the rate-of-deformation tensor D pl , which we 
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will subsequently call the plastic-strain rate as a matter of convention follows from the STZ 
dynamics. For a monotonically loaded, athermal system where there are no rate-dependent 
processes such as aging, which compete with the STZ-transition rates, and where we assume 
there to be a low STZ density, the flow rule can take the form, 

D pl = Fe -1/x , (1) 


where F = F(S) is a monotonic tensor-function of the deviatoric Cauchy stress S. One 
critical way the STZ theory differs from its predecessors (and indeed other theories) is 
through the introduction of the quantity \ and its relationship to an “effective temperature” 
T eS that is defined as 


Tpff — 


dU c 

dSr 


( 2 ) 


where U c and S c are the material’s potential energy and entropy respectively of only the 
configurational degrees-of-freedom, i.e. those degrees-of-freedom associated with the struc¬ 
ture of the material. This is to be distinguished from the usual thermalized temperature 
T which accounts for the degrees-of-freedom which relax on timescales short compared to 
the observation time. The typical definition of T is applicable to the fast, i.e. vibrational 
degrees-of-freedom, but the configurational degrees-of-freedom are typically out of equilib¬ 
rium. Some limited attempts to experimentally measure an effective temperature for disor¬ 
dered materials have been made [53|, as well as other direct, quantitative comparisons with 
experiments of bulk metallic glasses 48]. Atomistic simulations of two-dimensional glasses 
lave already suggested that T e g is linearly related to the local potential energy per atom 


54j, and earlier many arguments were made for the existence of the notion of an effective 


temperature, related to the entropy, that characterizes the number of phase-space paths 
available to the system |42|, l55l-l57]. 


In the effective-temperature STZ formalism, the dimensionless scalar field \ is defined as 
X — kBT e fi/E z , where ks is the Boltzmann factor. Although y is a dimensionless form of 
T cff we shall subsequently refer to it as simply “the effective temperature” for readability. 
Here E z is a typical energy required to create an STZ. In the athermal limit the dynamical 
equation for the effective temperature x takes the form 


. S : D pl . 

X =-(Xc 

^eff 


X) + D x V 2 x ■ 


( 3 ) 


The first term on the RHS in Eq. [3] represents a source of plastic work per unit time that does 
mechanical work on the structural degrees-of-freedom. The parameter c e g is the volumetric 
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effective-heat capacity with dimensions of energy per unit volume, determining the energy 
input per unit increment of effective temperature. In flowing regions y converges to a 
limiting value Xoo, which represents the steady-state effective temperature where the work 
done to shear the structure no longer causes an increase in disorder. In developing the 
analysis that follows we have assumed for simplicity that y^ is a constant for the entire 
system. In general Xoo is usually rate dependent [45[, however in the low-rate limit Xoo is 
approximately a constant. This assumption is consistent with the lower-limiting value of 
effective temperature found for simulations of sheared glasses in the steady state 


. The 


initial value Xo characterizes the structure of the gel in the pre-sheared state, and ideally 
would come from an analysis of the microstructural information of the system’s constituents. 
In the absence of this molecular-level information, the form of Xo including the mean value 
and any seeding of fluctuations about the mean are usually chosen in a way to best match 


the macroscopic behavior, e.g. the stress-strain curves of the material 


44 


46 


481 ]. Although 


the range of values of Xo for a particular system is significantly restricted by the nonlinear 
form of Eq. [3] and its stability. 

The final term in Eq. [3] describes the diffusion of the effective temperature through an 
effective thermal diffusivity D x = k e g/c e s with dimensions length-squared per unit time, 
where k e s is the effective thermal conductivity. Here we consider the simplified case where 
D x is a constant and diffusion occurs only in the presence of changes in the gradient of y. 
The diffusivity D x itself sets a lengthscale for the simulation that is approximately the size 
of an STZ. We plan to address the general case of the possibility of a strain-rate-dependent 
diffusivity in a subsequent paper. 

The effective-temperature theory as a rheological model is established in Sec. HH Sec¬ 
tion [Hi] presents the results of the small-gap effective-temperature model, making direct and 


quantitative comparisons with the recent experimental results of 


26J, including stress, flu¬ 


idization time, and velocity data. The unique two-state fluidization that occurs during the 
transient is described and explained in terms of the dynamics of the effective-temperature 
constitutive law. The robustness of the model, as parametrized for a particular applied 
shear rate, is discussed in terms of the power-law behavior found by the experimenters. We 
conclude in Sec lIVI with a discussion of the model and its implications for understanding 
transients during plastic flow in YSFs, and comment on the on-going development of this 
model. 
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II. BASIC THEORY 


In this section we adapt the effective-temperature definition taken from the STZ theory 
for a small-gap Couette-ccll system, whose geometry is described in Fig. QJi. The symmetry 
of the Couette cell allows us to treat the problem as an effectively one-dimensional system. 
We deliberately move to the small-gap limit where the curvilinearity vanishes as shown in 
Fig. Hb- This is consistent with experimental findings of a negligible inhomogeneity in the 


stress field across the gap such that the stress was reported as being uniform 


26 


29]. The 


complete absence of any stress gradient in the model allows us to test whether the transient 
shear banding can be accounted for purely by a microstructural heterogeneity occurring near 
the moving rotor. Indeed, we find that the effective-temperature theory is able to reproduce 
the experimental phenomenon in nearly every detail without including the stress gradient. 
We further assume that we are in the inertial limit, namely V • S = 0 which is also consistent 
with the experimental conditions. 

A key experimental observation of the rheology of simple YSFs is that the steady-state 
behavior of the system is described by the well-known Herschel-Bulkley (HB) relation 


26 


29], and this will be used to determine a flow rule for D pl . 


A. Plastic-Strain Rate 


The one-dimensional plastic-strain rate e pl from the STZ theory takes the form 


i pl = /Me- 1 '* . 


( 4 ) 


We next assume that the steady-state behavior of the gel is given by a HB form, namely 
s = s c Y A (e pl )" for some parameters n and A which characterize a system that begins to 
unjam and flow plastically when the deviatoric stress s reaches s c , the critical stress. Here 
the HB coefficient A has units of Pa-s n , and sets the timescale for the fluidization of the gel. 
The plastic-strain rate is related to the measure of disordering by a function / = /(s). We 
require that Eq. Q] a form of the constitutive law in the STZ theory, reduce to a HB relation 
when y = Xoo- Therefore e pl in the steady state must be 


cP 1 



( 5 ) 
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FIG. 1. a. The edge of the inner Couette-cell cylinder (i?i = 23.9 mm) and outer cylinder edge 
(R 0 ) with gel (white) where the gap width w = R 0 — R\ = 1.1 mm. The inner cylinder is held 
at constant velocity v a = 0.77 mm/s. The nucleation of a shear band begins at Ri and then 
broadens outward in the y-direction. b. The Couette-cell system in the small-gap limit where the 
curvilinearity vanishes, transforming the Couette-cell problem to a simple shear geometry. Here 
the system is fixed at y = w and driven with speed v a at the boundary y = 0. 


This allows / to be determined so that 

( \ i/ n 

and we find 

\ 1 /n 

■ ( 7 ) 

Upon substituting the expression for / into Eq. Q] we arrive at the expression for the plastic- 
strain rate, namely 

( \ l/n 

gi/xcc-i/x _ (8) 

This result provides a continuum constitutive-law relating the plastic-strain rate to the 
deviatoric stress within the framework of an effective temperature. 

B. Stress Rate 

To derive the differential equation for the deviatoric-stress rate we express the strain rate 
as the sum of the elastic and plastic components, namely 



e — e el + <i pl . 


( 9 ) 




















The velocity v = v(y) as a function along the Couette-ccll radius y in Fig. Q]is then given 
by 

v(y ) = [ e dy = [ e el dy + [ e pl dy ( 10 ) 


and when integrated across the full gap-width, one finds 


sw 


Vo = 


s — s c 

A 


l/n 


d/Xo 


rRi 


3 -!/x 


dy 


( 11 ) 


d \ ^ / JRo 

where y is the shear modulus (with dimensions of stress), and 7 = v Q /w is the average 
strain-rate across the gap found by imposing a velocity v Q at the inner cylinder. The shear 


experiment reported in 


26], to which we are seeking to compare in detail, uses a 7 = 0.7 


s -1 . We use this rate in all of our simulation results, except for the stress-time profiles of 
Fig [7] where each curve corresponds to a unique 7 . Because plastic flow is the dominant 
flow process for the transient shear band formation reported in [26|, we have assumed for 
simplicity that the elastic response is linear and that all non-linear behavior results from the 
plastic response. While the ramp-up of the elastic response is very brief, linear-elastic strains 
persist throughout the simulation and act to drive the plastic flow. The shear banding itself 
is an inherently plastic phenomenon and the relevant timescale for the transient dynamics 
comes from the plastic response. Rearranging we find that 


s — y 7 — p 


A 


l/n 


gl/Xoo_ 


"Ri 


3—i/x 


w JR 0 


dy . 


( 12 ) 


III. MODEL SIMULATIONS 


The equations of motion (EOM) of the effective-temperature model that we have pre¬ 
sented in Sec. UTlare 

/ \ l/n 

/ Q — Q \ ' 

(13) 


A 


g P l = ( ) gl/Xoog-l/x 


s = yi- y 


S — Sr 


A 

l/n 


l/n 


i r w 

gi/x°o_ / g—Vx 

W J 0 


X = — ( L ~A) e 1/x ~- 1,x (Xoo ~X) + D, 

Ceff V A J 


dy 

d 2 X 
x 8y 2 ‘ 


(14) 


(15) 


The three EOM take the form of coupled, nonlinear, partial integro-differential equations 
which we integrate in time. At time t = 0 , s = 0 so that the fluid is initially unstressed in 
our analysis. Boundary conditions for y must also be placed at the rotor and the outer wall 
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PARAMETERS 

SYMBOL 

UNIT 

VALUE 

Herschel-Bulkley exponent 

n 

- 

0.55 

Herschel-Bulkley yield stress 

Sc 

Pa 

29.4 

Herschel-Bulkley coefficient 

A 

Pa-s n 

1.60 

Diffusivity 

D x 

9 —1 

mm s 

10~ 6 

Shear modulus 

V 

Pa 

95 

Volumetric effective-heat capacity 

Ceff 

Pa 

5 xlO 3 

Steady state x 

Xoo 

- 

0.064 


TABLE I. The parameters characterizing the small-gap effective-temperature model’s constitutive 
law are presented with the values found to best match the stress-time and velocity measurements 
of the experiments. The denotes the parameter is dimensionless. A, n, and s c are from the 
Herschel-Bulkley rheology. The c e fj, D x , and Xoo are from the effective-temperature constitutive 
model. 


of the Couette cell. We impose a no-conduction boundary condition dx/dy = 0, since x is 
restricted to the gel itself. 

As discussed in the Introduction the initial value of the effective temperature y 0 should 
ideally follow from a detailed analysis of the microstructural information about the gel. 
Without such detail, we are left to conjecture an appropriate y D . We began by following the 
work of Manning et al. 44}, where we assume the initial condition can be decomposed into 
a background term and a fluctuation term about the background 


Xo = x“ G + $Xo 


( 16 ) 


To simulate the disordering of the gel near the inner cylinder (rotor) we chose the simplest 
type of fluctuation term—a perturbation near the inner cylinder. Equation [T6l which 
is shown in Fig. [2] was given the explicit form of a spatially uniform background y® G = 
0.041, \/y : 0 < y < w with a perturbation 5y 0 = 0.001, My : 0 < y < 0.03mm, near the 
inner cylinder where disordering of the structure would likely first occur, and thus nucleate 
a shear band [26]. We found that this value for y^ G generated a stress-time curve that 


reasonably matched the range from the experiment, and the value of 5y 0 was as small as 
possible without inducing purely homogeneous fluidization. The surfaces of both the inner 
and outer cylinders also present sites for structural heterogeneities. The slightly larger stress 
at the rotating inner cylinder favors localization there. Indeed, a boundary layer is known 
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to form at the inner cylinder in the experiments, although we do not explicitly model the 
layer here. As discussed earlier, we ignore the stress gradient in our numerical analysis as 
it is assumed to be negligible for a small-gap width, other than to preference the nucleation 
at the inner cylinder. The ability of this fluctuation at the inner rotor to grow and lead to 
strain localization in the form of a shear band depends on both the size of the perturbation 
and the average value of y D , in addition to the chosen constitutive law. It was shown in detail 
that small-amplitude perturbations can possibly grow during the transient regime if \o and 
its fluctuations meet an explicit criterion 44], However, even if they grow, implying that the 
equation of motion for x is unstable over some time interval, it does not guarantee that the 
fluctuations will form a shear band. According to 44], the magnitude of \ Q must be below a 
limiting value for linear stability during the transient regime in order for strain localization 
to be possible. Furthermore the \o chosen here has a numerical localization number which 
is within the range corresponding to at least partial to full strain-localization. We found this 
Xo field matches experimental width measurements of the shear band at early times, while 
also keeping the perturbation at the inner wall small enough to be a linear instability. Since 
larger values of the perturbation only enhance localization, we have deliberately chosen a 
value that would conceivably test the growth of a linear instability into a fully developed 
shear band. Alternatively a universal criterion for shear band formation in time-dependent 


fluids has been proposed by 


31] whereby the shape of the stress-time curve determines when 


shear banding can initiate. In particular the condition that <9 7 s + 7 d^s < 0 was calculated 
to hold Vf : t > 3 s. The first noticeable growth in x (but not yet reaching Xoo ) occurs at 5 
s and is shown in Fig. [2j exactly after the linearly elastic regime ends and flow is initiating 
in the form of a shear band. This criterion is therefore consistent with the onset of shear 
banding in the effective-temperature model. 


Table U contains the values of the parameters which appear in the small-gap model’s EOM. 
It includes those parameters which arise from the effective-temperature phenomenology of 
the model c e g and D x , as well as those associated with a HB Law A,n, and s c . Once the 
initial effective temperature field Xo was chosen and found to produce results which were 
within the range of the experimental stress-time profile, the parameters of the model were 
fit to the experimental velocity profiles and fluidization time, and to the precise inflection 
and shape of the stress-time curve. The central parameter in the effective-temperature 
hypothesis is the volumetric effective-heat capacity c e g- which appears in the equation for 
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FIG. 2. The effective temperature x at various times during the shear-induced deformation of the 
small-gap model simulation. The initial field x(2/? 0) = Xo (heavy black dashed line) containing a 
perturbation near the inner wall (y = 0) into a shear band as time evolves. The increasing value of 
X shows the band broadening across the gap (y-direction) before reaching a uniform steady-state. 
A strain rate 7 = 0.7 s -1 is applied at the rotor. 



FIG. 3. Plastic-strain rates at different times across the gap during the small-gap model simulation. 
As plastic strain develops in the gel, e pl takes on a large value near the shearing rotor (y = 0 ). As 
a shear band forms and traverses the gap the system fully fluidizes via a two-state process, until 
steady state is reached and everywhere e pl = 0.7, the value of the applied shear-rate. 

X■ The quantity c e g has the physical significance of being the amount of plastic work per 


12 































unit volume required to cause a fractional relaxation of x to its steady state, and its value 
was found to significantly determine the fluidization time and the abrupt transition of the 
entire system to the steady state. The shear modulus /i was set to match the experimental 
stress-time curve. 

The experimenters also reported a flow curve in a separate set of experiments from which 
they extracted values for a HB steady-state rheology. The values of these experimentally 
reported parameters, A = 9.8 (Pa s n ), n = 0.55, and s c = 26.9 (Pa) are inconsistent with 
the experimentally reported “steady-state” condition in the shear experiment that was used 
to produce the stress-time curve in Fig. [2 This is in part due to a waiting time of 100s 
before each stress measurement was made in the flow curve experiment, compared to the 
shear experiment, which we are modeling here, where the steady state does not emerge until 
10,000s. The values of A and s c are also likely dependent upon the particular sample and 
its preparation while the exponent n is believed to be more robust. Our constitutive law for 
the plastic-strain rate assumes some HB form must also hold for the longer-time protocol 
in the shear experiment, and hence cannot account for this discrepancy between the two 
types of experiments. As a starting point we initially used the HB value for A reported 
from the flow-curve experiments and then adjusted it to match the long-time behavior of 
the experiment plotted in Fig. [dj from which the value of s c was extracted. We also confirm 
the robustness of the exponent n, as the same value of n = 0.55 reported experimentally 
also best matched the data in the model. 

In the shear experiment where the rotor’s velocity was held at v = 0.77 mm/s, the gel was 
initially prepared to remove memory effects using the same protocol before each experiment 
26]: A pre-shear lasting 60 s at 10 3 s _1 was applied in the clockwise direction. Then 60 s at 
10 3 s^ 1 was applied in the counterclockwise direction. The shear was then instantaneously 
stopped and the gel was allowed to rebuild for 120 s. The experiments then involved shearing 
at a constant rate for 10 4 s while the stress was measured. A shear band was observed to 
nucleate near the rotating inner wall of the cell. After nucleation the band widened until the 
onset of fluidization, at which point the entire gel transitioned to a homogeneously flowing 
state, revealing a distinct process of unjamming. During the transient characterized by the 
broadening of the shear band, the shear rate and velocity profile of the gel outside the band 
were found to be non-zero. 

The gel’s restructuring under shear is captured by the evolution of x during the numerical 
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FIG. 4. The deviatoric shear stress as a function of time for the small-gap model and the experi- 

261 ] under an applied strain-rate 7 = 0.7 s~ x at the rotor. Complete 


mental stress measurements 


fluidization of the gel is experimentally observed between 2500-3000s, and is well captured by the 
simulation. 


simulations. In Fig. [2] y is plotted at different times during the deformation. The bottom¬ 
most, dashed black curve in Fig. [2] illustrates x(y, 0) = Xo, the state before shearing begins. 
The time-evolution of the plastic-strain rate is described in Fig. [3] and reflects these changes 
in x- Central to the theory is the notion that the local plastic-strain rate depends on the 
changing state of the amorphous structure, which is described by x — x{Vi t) ■ At any time, 
X quantifies the disordering of the material as it is strained or stirred. In the athermal limit, 
the time-evolution y will be nonzero only if a state of stress greater than a critical stress 
is present or if there are spatial variations large enough to activate diffusion through V 2 y. 
Stresses above the critical stress induce microstructural rearrangements which cause x to 
evolve towards its steady-state value Xoo- 

In our effective-temperature model the steady state is determined by the two terms on 
the RHS of Eq. [31 a source term S : D pl (y^ — x) / c eff and a Laplacian term -D x V 2 y which 
acts to diminish gradients in y. As such, the steady state simply occurs when y = y^, since 
the Laplacian of y becomes negligible in the steady state. 

As was observed experimentally, the simulations reveal the formation of a shear band, 
which nucleates from the perturbation near the inner cylinder and begins to grow, broadening 
across the gap. I 11 the simulations this is followed by an abrupt transition during which 
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the gel reaches a fully unjammed state and completely fluidizes. We found the value of 
Xoo = 0.064 to best match the sudden transition to homogeneous flow. The stress-time 
curve of the simulation also agrees with the experimental measurements as seen in Fig. 01 
where in the model the gel undergoes an extremely brief, linearly elastic regime followed 
by plastic flow (including the transient shear banding) until the steady-state stress given 
by the HB relation is reached. The time needed for the shear band to grow sufficiently to 
accommodate the strain rate imposed across the gel results in the overshoot in the stress 
that is seen in Fig [4] This is related to the value of the initial condition for the effective 
temperature. For example, a stiffer, stronger, more ordered system (lower y D ), will have a 
larger stress-overshoot. 

The sudden fluidization of the gel after the formation of the shear band is also evident in 
the velocity profiles shown in Fig. [5] The portion of the velocity curves with the significantly 
steeper slope near the inner cylinder is indicative of a shear band, while the lower slope near 
the outer cylinder suggests that the gel outside the shear band is also gradually flowing as 
it fluidizes. Experimental measurements confirm the non-zero velocity for the gel in front of 
the shear band. Direct, quantitative comparisons of the velocity profiles of the simulations 
and the published experimental data are difficult because of the significant wall slip found 
in the experiment especially at early times, effects which are not included in the model. In 
Fig. [5] we have compared three instances from the simulations against the experimental data 
with reasonable agreement as the fluidization time is approached, experimentally between 
2500-3000 s. 

We see the presence of the two distinct fluidization processes in the model as a re¬ 
sult of the effective-temperature dynamics and the particular form of the constitutive law 
from the STZ theory that we are using. The equation of motion for y takes the form 
of a second-order parabolic partial differential equation (e.g. the heat equation) with 
source and Laplacian terms as seen in Eq. [3] or [151 In the case of spatially homoge¬ 
neous solutions, like those found in the steady state, Eq. [15] reduces to y — 9x(x)> where 
9\{x) — ~ e 1 /* 00-1 ^ (Xoo — x) is a nonlinear, inhomogeneous term in y. Here, 

the function g x {x ) is similar in form (although asymptotically different) to the logistic re¬ 
action term found in the F-KPP (Fisher-Kolmogorov-Petrovsky-Piscounov) equation that 


describes solidification and reaction-diffusion fronts 


59|-6l|. In the case of the F-KPP equa¬ 


tion, a consideration of the stability of equilibria reveals that there is one unstable fixed-point 
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FIG. 5. Velocity profiles at different times during the small-gap model simulation (solid curves) 

compared to experimental measurements at times 27 s (— • —), 1,730 s (-), and 2,250 s (•■•)■ 

The velocity evolves from a discontinuous profile over the gap, indicating the coexistence of a shear 
band with more gradual fluidization in the material outside the band. There is significant wall slip 
at early times in the experiments. A nearly linear profile in the steady state is observed in both 
the experimental measurements and the simulations. 


corresponding to a metastable state, and two stable fixed points into which material can 
be transforming as illustrated in Fig. | 6 ] Similarly, g x also has a stable fixed-point which 
corresponds to the shear-banding region where y has reached the steady-state value Xoo- 
But instead of a metastable state outside the shear band, g x —> 0 exponentially, as y —>■ 0. 
Consequently the gel outside the shear band is extremely sluggish, although it remains sub¬ 
ject to gradual fluidization in the presence of a shear stress, even far from the shear-banding 
region. Ultimately it is the competition between the slow fluidization of y outside the shear 
band and the forward growth of the shear band that begins at the inner cylinder, that acts 
as the crucial element in the effective-temperature model leading to the sudden transition 
to homogeneous flow. 

From Fig. [I] we see the simulation reaches the steady-state stress after approximately 
2,600 s. This is also observed from the velocity profile across the gap shown in Fig. 0 The 
experimentally reported fluidization time for an applied shear-rate of 7 = 0.7 s -1 was also 
approximately between 2500-3000 s. The experiments report a power law for the fluidization 
time Tf ~ 7 -2 - 3 as a function of the applied shear rate 7 over a range 0.08 < 7 < 10 (s _1 ) 
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FIG. 6. Fluidization outside the transient shear band: A comparative sketch of the form of the 
F-KPP equation dtu — d%u = g u [u) (dashed curve) and d t x ~ D x d%x = 9\(x) (solid curve) in 
the effective-temperature theory. Points 1 and 3 represent stable fixed-points for g u (u ) which 
can describe distinct stable phases of a system into which material that is in a metastable state, 
represented by the unstable fixed-point at 2, could be transforming. In the effective-temperature 
model however, only one fixed point exists at 4; it is stable and corresponds to the shear band 
where x = Xoo> the flowing steady state. For small x near point 1, corresponding to the gel outside 
the shear band, g x {x ) —>• 0 due to e -1 /*. However no fixed point exits at 1 for <? x (x), and this 
feature allows the gel outside the shear band to slowly fluidize through a distinct process of its 
own. 

using several boundary conditions and sizes of gap width. The effective-temperature model 
is able to recover approximately the same exponent for a subset of the rates from about 
0.1 < 7 < 2.0 centered on 7 = 0.7, the specific 7 for which the model’s parameters have 
been determined with w = 1 . 1 mm and where a no-slip boundary condition is implicitly 
assumed. The faster the shear rate in the model is, the faster fluidization and homogeneous 
flow is reached, and consequently the transient shear band travels a shorter distance across 
the gap. Figure 0 shows the power-law fit for our simulations. Outside this subset the 77 
predicted by the model begins to deviate by about an order of magnitude. Interestingly 
however, the experimental fluidization times outside this subset can still be recovered if 
c e ff is adjusted (re-parametrized) to a lower value for very small shear rates, and increased 
for the largest shear rates. This would imply that the fraction of plastic work that goes 
into disordering the gel increases (decreases) as the strain rate increases (decreases), and 
the power law of the fluidization time 77 ~ j~ a may arise in part from this aspect of the 
physics. This could reflect rate effects in the gel’s mechanical response which cause the 
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FIG. 7. a. The simulated fluidization times ry of the model which was parameterized for the 
case 7 = 0.7 s -1 where w = 1.1mm, and assuming a rough boundary condition. Experimenters 
have reported a universal power law ry ~ y -2 - 3 for a wider range of shear rates, various boundary 
conditions, and several different gap widths, b. The stress-time curves at different imposed strain- 
rates (inner cylinder velocities) from which the fluidization times (at the onset of the steady state) 
were extracted. 


gel to accommodate deformation through structural disordering at high rates that it could 
accommodate through independent relaxation processes at low rates. An analysis of the 
effect of strain-rate dependence is beyond the scope of the theme of this paper, since the 
model here has not explicit rate-independent. We nonetheless think this is an important 
direction for future work. To our knowledge no other theoretical models describing the 
solid-fluid transition make any quantitative predictions for the fluidization times which agree 
(within an order of magnitude) with these experiments over any range of rates. 
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IV. CONCLUSIONS 


We have presented a phenomenological effective-temperature model for the small-gap 
limit of Couette-cell shear experiments of a carbopol gel where transient shear banding 
is observed. This theoretical model is based on the STZ description, and results in two 
distinct fluidization processes: the shear band formation and the simultaneously competing 
homogeneous fluidization. We have made direct, quantitative comparisons with experiments 
which demonstrate reasonable agreement with the stress-time behavior and velocity profiles 
of the gel. 


In this theory, consistent with the interpretation of 


26], stress gradients are not the 


primary cause of the strain localization or the fluidization. The small-gap limit of the 
effective-temperature model strongly suggests that the transient shear banding and sudden 
fluidization is primarily a result of microstructural disordering originating from structural 
heterogeneities at the inner cylinder of the Couette cell where the rotor applies the shear. 
While the gradient is at least negligibly present kinematically, our analysis suggests that it 
may play little role in the constitutive relation for a YSF under shear—and no role for a 
small-gap system. This is in contrast to other theoretical descriptions in which the strain 
localization is directly associated with the fluidization and comes solely from the non-linear 
response of the fluid coupled to the stress inhomogeneity that arises from the cylindrical 


geometry of the Couette-cell 


32], One possible experimental check woulc 


be to see if such 


26f were repeated 


transient shear banding and fluidization would take place if the experiment 
for a different shear geometry, such as shearing the gel between two parallel plates. The 
absence of any normal stresses from the model’s constitutive law suggests to us that such 
stresses are not necessary to account for the transient shear banding, as was also found by 
the A-based rheological models discussed in the Introduction. 


The effective-temperature theory also improves several aspects of existing visco-plastic 
models which lack any notion of a yield stress, a defining physical feature of YSFs. The 
mechanical response of the gel is quantitatively well captured by the small-gap simula¬ 
tions, as evidenced by the stress-time curve. The small-gap simulations in this paper reveal 
the two-stage fluidization seen in experiments, and we have shown that in the model it 
is directly attributable to the particular form of the constitutive law postulated by the 
effective-temperature theory. The parameters of the theory provide a physical connection 
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to the thermodynamics of the gel’s structural state under shear. One open question is the 
origin of the strain-rate dependence of the fluidization time for very large and very small 
rates . A possible explanation is that the volumetric effective-heat capacity, describing 
the fraction of plastic work that disorders the gel and emerges phenomenologically in the 
effective-temperature dynamics, is strain-rate dependent and can no longer be approximated 
as uniform over a wide-enough range of rates. 

The effective-temperature approach that we have proposed offers an important rheological 
model for both current and future experiments, and responds to a need that was identified 
by the exp erimenters for developed theoretical models of transient shear banding in simple 
YSFs [26]. We also believe our approach would benefit from further development, and 
are currently investigating an extension of this approach for finite-gap systems with non- 
negligible stress gradients. We are also currently attempting a version of this theory with 
a rate-dependent diffusivity and steady-state effective temperature. Beyond this a number 
of consequential phenomena exist which we have largely ignored in the present model, e.g. 
aging and wall slip, both of which are known to play at least some role in the recent shear 
experiments of carbopol. A generalization of this model and incorporation of such additional 
physics is also a matter of our on-going research. 
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